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Abstract 

Minimally doubled fermions have been proposed as a cost-effective realization of chiral sym- 
metry at non-zero lattice spacing. Using lattice perturbation theory at one loop, we study 
their renormalization properties. Specifically, we investigate the consequences of the breaking 
of hyper-cubic symmetry, which is a typical feature of this class of fermionic discretizations. 
Our results for the quark self-energy indicate that the four-momentum undergoes a renor- 
Qh| malization which is linearly divergent. We also compute renormalization factors for quark 

bilinears, construct the conserved vector and axial- vector currents and verify that at one loop 
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the renormalization factors of the latter are equal to one. 
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1. The past ten years have witnessed two major breakthroughs in lattice QCD. The first con- 
cerns the significant acceleration of simulation algorithms for dynamical quarks, in particular 
as the light quark masses are tuned towards the chiral regime. Owing to these developments, 
simulations with pion masses close to the physical value have become routine. The second 
achievement was the solution to the long-standing problem of constructing discretizations 
of the quark action which preserve chiral symmetry, and the realization of the role of the 
Ginsparg- Wilson relation [1,2]. However, the well-known discretizations based on the perfect 
action formalism [3], or, alternatively, the domain wall [4,5] and overlap constructions [6] all 
involve non-local interactions. As a result, their implementation is numerically much more 
expensive than conventional Wilson [7] or staggered [8] fermions. 

Minimally doubled fermions [9-13] share the desirable features of strict locality with tra- 
ditional discretizations, whilst preserving exact chiral symmetry for a degenerate doublet of 
quark fields. The question whether or not they are suitable for the determination of hadronic 
properties in practical simulations has not been thoroughly investigated so far. Before em- 
barking on extensive numerical studies of minimally doubled fermions, it is useful to examine 
some of their properties in perturbation theory. In this letter, we present our results for 
the quark self-energy and the renormalization properties of quark bilinears at one loop in 
lattice perturbation theory. In particular, we shall elucidate the consequences of the break- 
ing of hyper-cubic symmetry, which is a typical feature of this class of lattice actions. Our 
main findings indicate that hyper-cubic symmetry breaking generates a renormalization of 
the quark's four-momentum. 

After fixing our notation in section 2, we list expression for the quark propagator and the 
vertices in section 3. The perturbative calculation of the quark self-energy at one loop in lattice 
perturbation theory is presented in section 4. Sections 5 and 6 discuss the renormalization 
properties of quark bilinears and provide expressions for the conserved vector and axial- vector 
currents. In section 7 we present our conclusions and discuss the consequences of our results 
for numerical simulations. 



2. Following the works of Borigi [12] and Creutz [13] we employ the particular construction 
of minimally doubled fermions of ref. [13]. The corresponding lattice Dirac operator respects 
chiral symmetry, and is also 0(a)-improved. For massless quarks, the general expression reads 

V = D + D -2iT. (1) 

In momentum space the terms D and D are given by 

Dip) = i 5^(7^ sinpf,), D(p) = i cosp^). (2) 

n n 

The matrices T and 7^ are defined by 

n 

with r 2 = 1. Other useful relations are 

E^ = Ei = 2r < { r >>}=i> {r,v„} = i. (4) 

The construction of the Dirac operator T> involves a particular linear combination of two 
(physically equivalent) naive fermion actions, corresponding to D and D ineq. ([I]). The first 
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term, as is widely known, has 16 zeros in the first Brillouin zone, when any component of p is 
equal to or n. The term D has also 16 zeros, which are, however, positioned at the points 
where p^ = ±7r/2. 

As was shown in ref. [13], the presence of the term —2iT in eq. ([1]) guarantees that T>(p) 
exhibits only two Fermi points, located at p = (0,0,0,0) and p = (vr/2, tt/2, tt/2, vr/2), which 
represent two degenerate fermion species of opposite chirality. Here we simply note that the 
extra zeros of D at the corners of the Brillouin zone are lifted by the presence of D — 2iT. An 
entirely similar statement applies to the zeros of D. 

The matrix T is not unique: There are altogether 16 definitions of T, in which the 
coefficient of a particular gamma-matrix can be chosen as 1 or —1. Each choice selects 
different zeros of P, which correspond to the physical degrees of freedom, but otherwise all 
such definitions yield an equivalent theory. 

The inclusion of the term proportional to T implies that the action is no longer symmetric 
under the full hyper-cubic group. Depending on its definition, the matrix T selects a particular 
direction in Euclidean space. The action of minimally doubled fermions is only symmetric 
with respect to a subgroup of the full hyper-cubic group, which preserves this fixed direction 
(up to a sign). For the action considered in this paper, which corresponds to the definition of 
r in eq. ([3]), this is the positive major diagonal. For other minimally doubled actions, such as 
those considered in refs. [9,10], it is the temporal axis. 

The breaking of hyper-cubic symmetry implies the possibility of mixing with operators 
of different dimensionality, such as tpTOip. In refs. [14-16] it was argued that mixing with 
dimension-3 operators cannot be avoided. In any case, the lack of hyper-cubic symmetry gen- 
erates not only mixing with dimension-3 operators, but also mixing with marginal operators 
of dimension 4. In the next two sections we will show that hyper-cubic symmetry breaking 
generates a linearly divergent additive renormalization of the quark's momentum. 



3. By inverting the Dirac operator of eq. ([I]), and restoring the proper factors of a, we obtain 
the fermion propagator 

— i 7At(sin ap^ - cos ap^) - «r(^J cos ap^ - 2) + am 

S(p)=a— ^ 1 -. (5) 

y sin ap^ cos ap v — 2 sin ap^(cos ap fl + 1) — 2 cos ap^ + 8 + (omo) 

H v 

Unlike many standard fermionic discretizations one finds that the denominator of this prop- 
agator cannot be cast into a form which possesses a definite behaviour under parity transfor- 
mation in each single coordinate (pi — > —pi). This is not surprising in view of the fact that 
the definition of V singles out an intrinsic Euclidean direction. 

By using the identities {7^,7^} = {7^>7^} = 2 <W and {7m>7^} = 1 ~ 2 <W> tne above 
propagator can also be written in a form which is more convenient for lattice perturbation 
theory, i.e. 



^ [7/i sin apf, - 2 7^ sin 2 ap fJi /2^ + am 



S(p) = a = ^1 > (6) 

4^ sin 2 ap^/2 + sin apA sin 2 ap^/2 - \ ^sin 2 ap u /2j + (am ) 2 

where the limit of small p (i.e. the continuum limit) becomes more transparent. In particular, 
one can use the standard methods for the treatment of lattice divergences, and calculate for 
example the so-called J and I — J parts along the lines of [17]. 
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Via the substitution ap^ — » ir/2 + ap^ one obtains the propagator for the fermionic mode 
associated with the Fermi point at ap^ = ir/2: 



-i 



^2 [ - sin ap^ - 2 7 M sin 2 ap M /2j + am Q 



S'(p) = a = 1 ^ . (7) 

4 sm2 a P^/^ ~ sm a Pfj. ( sm2 — \ sm2 a Vvl^ ) + (a^Tio) 2 



By changing the direction of the four- momentum p^ and exchanging 7^ with 7^ , one recovers 
the propagator of eq. (J6]). Since 75 = —75 this implies that the modes corresponding to the 
two Fermi points have indeed opposite chirality. 
The quark-quark-gluon vertex is derived as 

r r 1 \ • ( a(p 1 +p 2 ) fl , . a(p l +p 2 )A 

Vi{pi,p 2 ) = -igo I 7m cos j ~7 M sm J , (8) 

and the quark-quark-gluon-gluon vertex is 

V 2 (pi,P2) = 7; ia 9o ilfn sm ^ + 7 ^ COS 2 J ' ^ 

where pi and j>2 are the incoming and outgoing quark momenta at the vertex. Following 
ref. [18] , the expressions for the vertices can be easily derived by comparing the Dirac operator 
of minimally doubled fermions, eq. (pQ), with the Wilson-Dirac operator 

D w {p) = - |i7 M sinap M + r (1 - cos ap^)j + m , (10) 

and noting that the hopping terms of these two actions are related by the replacement r — ► 
— ij'. Indeed, since the terms e iaPti and e~ taPti are coupled to the Fourier transforms of U^(x) 
and Ujt(x — ajx) respectively, it is sufficient to substitute r — > —27^, in order to obtain the 
vertices for minimally doubled fermions from those of the Wilson case. 



4. We are now going to describe the calculation of the quark self-energy at one loop. FigureQ] 
lists all diagrams which are relevant for the perturbative calculation in this letter. 

Using the expression for the vertex V 2 (p,p), the tadpole contribution to the self energy 
(diagrams (g) and (h) in Fig. [1]) is easily computed. In a general covariant gauge, where 
d^A^ = 0j the expression is 

h ' t i 1 ~ \ {l ~ a) ) ' iag2 ° CF S + ( r - + °( a2 )) 

= ^^(l-i(l-a)) (^ + i^(r- 7M ))+0(o), (11) 

where Zo is given by [19-21] 

/""/a ^4„ 1 14 /an \ 

., /o (WF = al549333 ''- = 2 " 66100 i^' ?2 = ;?I> 2 (ir)' < 12 » 

Terms of 0(a) and higher are not important here. Since 7^ = 2r, the result of the one-loop 
tadpole is 

9lc F ^ (1 - |(i -«))(* +^)- (13) 
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(a) vertex (b) sail 




(c) sail (d) operator tadpole 




(e) leg self-energy (f) leg self-energy 

(excluding tadpoles) (excluding tadpoles) 




(g) leg tadpole (h) leg tadpole 



Figure 1: The diagrams needed for the one-loop renormalization of the lattice operators. 



The term proportional to ifl is the same as for Wilson fermions, while the other term, as 
already noted in [14] , would imply a power-divergent mixing of order 1/a with the dimension- 
three operator ipTi/i, provided that there is no cancellation by an analogous term coming from 
the contribution of the sunset diagrams to the self-energy (diagrams (e) and (f) in Fig. [TJ. In 
this section we show that there is no such compensation. 

We have computed the sunset diagram using special computer codes written in FORM 
[22,23] and Mathematica, and also checked it against calculations by hand. The result of this 
diagram is 



^sunset 



(p,m ) = n 



+ m --^C F 



16vr 2 



16vr 2 



C F 



log a V - 5.42642 + (1 - a) - log a V + 7.850272 



4 log aV - 29.48729 + (1 - a) ( - log a 2 p 2 + 5.7920K)) 



a 2 a 2 T 

+ 1.52766 • ArCfirVfc + (5.07558 + 6.11653(1 -a)) .Jzp-C F -i-. 

16tt z ^— ' I67H a 



(14) 
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Note that gauge invariance forces the terms proportional to (1 — a) to be the same as, for ex- 
ample, in the case of Wilson or overlap fermions. This is an important check of the correctness 
of our calculations. 

The total contribution of the one-loop diagrams to the quark self-energy is then 

S(p,m ) =iySi(p) + m £ 2 (p) +ci(gjj) ■ ^r^p^ + c 2 (go) (15) 



where 



=2(p) = Ac> 



log aV + 6.80663 



+ 



(1 - a) ( - log aV + 4.792010) 



(16) 



16vr 2 



4 logaV - 29.48729 



+ (1 - a) ( - log aV + 5.792010) 



(17) 



Cl (<7 2 ) = 1.52766 -^C F (18) 

c 2 ( ff0 2 ) = 29.54170 • jjL C F . (19) 

As indicated above, the two terms proportional to T/a arising from the tadpole and the sunset 
diagrams do not cancel — they actually reinforce each other. Note, however, that the parts 
proportional to (1 — a) cancel exactly, as required by gauge invariance. 
The full inverse propagator at one loop can be written as 

V m ) = (l - E X - |) • {<y + m (l - £ 2 + E x + |) - ^ £ W - ^ ij, (20) 

where we have collected all terms proportional to itf in the wave-function renormalization, 
which then contains, in addition to the standard term £1, also c\. By contrast, the linear 
divergence (i.e. the term proportional to c 2 ) must be absorbed into a redefinition of the 
four-momentum, which amounts to a uniform additive shift, 

PI) 

After replacing p by p' and neglecting terms of 0(<7g), we obtain 

*J» + Z m m - 1 5 ci (gg) 7/^ 
where the wave-function renormalization at one loop is given by 

*=(l-=i-!)"\ (23) 

while 

- (£ 2 - Sx - |) (24) 
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is the result for the quark mass renormalization factor. 

Our results demonstrate that the self-energy generates a power-divergent mixing with 
an operator of the form T/a. As can be seen from the Dirac structure, this mixing is not 
a renormalization of the mass. Indeed, chiral symmetry protects the quark mass against an 
additive renormalization like in the Wilson case. Rather, the power-divergent mixing implies 
that all components of the four-momentum p^ are shifted under renormalization by an equal 
amount, i.e. p M — > p'^ = p^ + const. /a. The constant can be determined either order by order 
in perturbation theory, or at the non-perturbative level in a Monte Carlo simulation. For 
instance, in our perturbative calculation it is given by — C2(<7q)/2- 

It is important to realize that the term proportional to ci(<7q) cannot be absorbed into a 
redefinition of p„, since otherwise the conserved vector and axial- vector currents do not have 
unit normalization. We address this issue in more detail in section 6. Thus, the renormalized 
quark propagator, eq. (f22j) . contains also a term X)„^„7/xJV This should not come as a 
surprise, since the present formalism is no longer isotropic. Note that the presence of this 
term does not move the pole of the propagator at p = 0. 

Since the mass is protected from an additive renormalization, the redefinition of the 
four-momentum amounts to a renormalization of the velocity. Noting that p^ is a fermionic 
momentum (it is the external momentum of the fermionic self-energy), and that no such 
phenomenon can occur for the gluonic self-energy, we can interpret this mixing as a renor- 
malization of the quark velocity. These findings support Creutz's conjecture [11, 13] that 
"interactions at finite lattice spacing can result in the gluons and fermions not having the 
same speed of light." 

We remark that the power-divergent mixing proportional to C2, as well as the one of the 
same dimensionality (proportional to ci), occur among operators which are not invariant under 
the hyper-cubic group. Such mixings are lattice artefacts which are peculiar to minimally 
doubled fermions. 

Once the above subtraction has been made, and p^ replaced with p'^ in all renormalized 
quantities, the power divergence disappears. Although we have no proof, it is not unreasonable 
to expect that this can be done consistently at every order of perturbation theory, similar to the 
subtraction of the 1/a divergence in the self-energy of Wilson fermions, which is consistently 
removed from the theory by the replacement too — > m q = — m cr . 

One may wonder how the redefinition of the four-momentum affects numerical simulations 
and how it could be determined non-perturbatively. We postpone this discussion to our 
conclusions in section 7. 

5. We have also computed the renormalization factors of local bilinears, and we list the 
results for the vertex diagrams below. The complete renormalization factors are obtained 
after including wave- function renormalization, which is achieved by adding the contributions 
of Ti\ and ci, eqs. (|16H and (]15p . of the self-energy. For the scalar density the result for the 
vertex (diagram (a) in Fig. [I]) is 



Here there is no mixing term arising from the breaking of hyper-cubic invariance. The only 
such mixing occurs after adding the wave-function contribution, which includes the term 




C F 



4 logaV + 29.48729 + (1 - a)( logo 2 p 2 - 5.792010 



') 



(25) 
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proportional to ci(<?q). For the vector current the vertex diagram yields 



9 ° r n, 



log aV + 9.54612 + (1 - a) ( log a 2 p 2 - 4.79201o) 



16vr 2 

where the coefficient of the mixing is given by 



+ c ; tx ( 5o 2 )-r, (26) 
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c v v tx (g 2 ) =-0.10037-^0. (27) 

This is a mixing with an operator of the same dimension, which is not invariant under the 
hyper-cubic group. Note that there can be no power-divergent mixing here (and in all the 
other bilinears), as one can see by simple dimensional counting. 

As a consequence of chiral symmetry, the vertex corrections are identical for the vector 
and axial-vector currents, and the same is true also for the scalar and pseudo-scalar densities. 
We have verified this in the course of our calculations. After taking the renormalization of 
the wave- function into account, the renormalization factors Zy and of the local vector 
and axial- vector currents are not equal to one. In order to identify the conserved currents, 
which are protected against renormalization, one has to consider the chiral Ward identities. 
We postpone this discussion to the next section. 

Finally, for the tensor current we obtain the result for the vertex diagram as 



9o n _ 



2.16548 + (1 - a) logoV - 3.792010 



(28) 



Again, the breaking of hyper-cubic invariance does not generate any extra mixing, apart from 
the one arising from the self-energy. 

We end this section with a brief comment on the renormalization of the quark mass. Chi- 
ral symmetry protects the bare quark mass tuq from undergoing an additive renormalization. 
The relation between the bare and renormalized quark masses, thq and tor, respectively, is 
then obtained via 

m R = Z m m , (29) 

where Z m is given in eq. (|24p . The full expression for the renormalization factors of the scalar 
and pseudo-scalar densities in perturbation theory at one loop is 

Z S = Z P = 1- (A s + E 1 + £L), (30) 

where the results for the self-energy contributions Ei, and c\ are given in eqs. (|16p and (|18[) . 
Here As is the result for the one-loop vertex diagram of the scalar density, given in eq. (|25[) . 
which is exactly equal to the 0(c/g)-contribution to the quark self-energy S2, but comes with 
an opposite sign. Thus, when we compare with eq. (|24p . we see that the renormalization 
factors Z$ and Zp of the scalar and pseudo-scalar densities satisfy 

1/Z m = Z S = Z P , (31) 

where the last equality is a consequence of chiral symmetry. We have thus verified at one loop 
in perturbation theory, that the renormalization of the quark mass for minimally doubled 
fermions has the same form as, say, in the case of overlap fermions. 
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6. One can derive expressions for the conserved vector and axial-vector currents via the chiral 
Ward identities along the lines of ref . [24] . The lattice action of minimally doubled fermions 
in position space reads 

Sf = °" E ^ E [^O^ + Wy)Up{x)1>{x + Oft) 

- ^(x + ajl)^ - i^)Ul(x)ip(x) +ip(x)(mo-^-^'ip(x) .(32) 



It is important to recall that the Dirac spinor ip in this expression describes a degenerate dou- 
blet of quarks. The action in the massless case is invariant under an axial U(l) transformation, 
and it is then clear that the chiral Ward identities associated with this exact symmetry yield 
the isospin-singlet currents of the theory. 

If one applies the usual vector and axial transformations, i.e. 



Syip = iayip, Syip = —itpay, 



one identifies the conserved vector current as 



V£( x ) = \ (^( x ) (7/* + U ^ x ) ^ x + a ^> + ^ x + (7m -il'n) U^x) ip(x) 
while the axial-vector current (which is conserved in the massless case) is given by 



A ^( x ) = \ H>(x) (7 M + 75 U^x) tfj(x + ap) + ip(x + a/2) (7 M - i^) 75 E^(x) if>(x) 



(33) 



(34) 



(35) 



Below we list the results for the individual diagrams of the conserved vector current. Due to 
chiral symmetry, the corresponding expressions for the conserved axial current are trivially 
obtained by replacing 7^ with 7^75, and T with 1^75. The vertex (diagram (a) in Fig. [Q) gives 
the result 



9l 



log aV + 0.61800 + (1 - a) ( log ay - 1.73375 



1671-2 

where the mixing coefficient c^* x is given by 



+ c?*(gi)-T, (36) 



-0.43749 • 



16vr 2 



C F . 



The result for the sails (diagrams (b) and (c) in Fig.[T]) is 



167T 2 



Cf7m 



4.80841 -6.11653(1 - a) 



+ C(ffo)-r, 



where <4 C S is obtained as 



-1.09017 



16vr 2 



C F . 



(37) 



(38) 



(39) 



Finally, the operator tadpole (diagram (d) in Fig.Q]) gives the same result as for Wilson 
fermions: 

Jin. - - I, { \ ( 40 ) 
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Summing up all contributions gives 



5 o r n, 



logaV - 6.80664 + (1 - a)( loga 2 p 2 - 4.79202) 



16vr 2 

where the total mixing coefficient is given by 



+ c vc (g 2 )-T, (41) 



CvcQ/o) = -1-52766 • -f^ C F . (42) 



These numbers exactly compensate the contributions of and c\ of eqs. (fl~6j) and (fl~8|) 

of the quark self-energy. Although we are not yet able to give an algebraic proof that the 
term i\c\ Y^u+v luPv i n the denominator of the self-energy, eq. ([22]) . cancels the unwanted 
contribution of the conserved currents proportional to T, we know that this has to happen, 
for otherwise the renormalization factors of the conserved currents would be different from 
one. We have explicitly derived these currents using chiral Ward identities, corresponding 
to transformations which leave the Lagrangian invariant, and this proves that these currents 
are conserved. Thus, it is certain that such cancellation must occur, and the one-loop result 
c vc = — ci (which holds to all significant digits that we have achieved) lends support to this 
statement. 

It is for the above reasons why the term proportional to c\ cannot be absorbed into 
the redefinition of the four-momentum, as this would spoil the correct renormalization of the 
conserved currents. Our one-loop calculation has thus confirmed within our numerical preci- 
sion that the renormalization factors of these currents is unity, as expected. It is remarkable 
that the use of the conserved currents exactly cancels not only the self-energy terms that 
contribute to the multiplicative renormalization, but also the mixing with dimension-four 
operators, coming from hyper-cubic symmetry breaking. Of course, radiative corrections to 
quark bilinears cannot generate any terms proportional to 1/a, and so an uncancelled, power- 
divergent factor arising from the self-energy remains in the total renormalization constant. 
However, as already stated, the latter can be absorbed into a redefinition of the renormalized 
four-momentum. By contrast, for the local vector and axial- vector currents any mixing com- 
ing from the hyper-cubic breaking remains uncancelled after the self-energy has been added 
to the vertex diagram. 



7. In this article we have presented the first perturbative study of a particular realization 
of minimally doubled fermions at one-loop order. Our analysis has shown that minimally 
doubled fermions are described at one loop by a fully consistent quantum field theory. We 
have elucidated the consequences of the breaking of hyper-cubic symmetry. In particular, we 
found that in the Borigi-Creutz construction [12,13], all components of the four-momentum 
undergo a subtraction under renormalization, which consists of a uniform shift. 
Furthermore, local vector and axial-vector currents can mix with other operators of the 
same dimension which are not invariant under the hyper-cubic group. By contrast, no such 
mixing occurs for the scalar density and the tensor current. We have derived expressions for 
the conserved isospin-singlet vector and axial currents, which involve only nearest-neighbour 
points. They do not undergo any mixing, and we have verified that their renormalization 
constants are one. In fact, apart from the staggered formulation, minimally doubled fermions 
are the only discretization which yields a simple expression for a conserved (point-split) axial- 
vector current. 

It remains to discuss the implications of our findings for practical simulations. Since 
there is one exact chiral U(l) (8>U(1) symmetry, there must be exactly one Goldstone boson as 
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the quark mass is tuned to zero. It is natural to associate the neutral pion with this particle. 
The charged pions, by contrast, will retain a finite mass in the chiral limit at non-zero lattice 
spacing. An interesting observation is that the influence of disconnected diagrams, which are 
required for the determination of the 7r° mass, must become weaker as the lattice spacing 
goes to zero, since the exact masslessness of the charged pions is recovered in the continuum 
limit 

Even though the neutral pion in the chiral limit is exactly massless in this theory, the 
renormalization of the four-momentum will modify the rate of the exponential fall-off of its 
two-point correlation function. Then, the extracted energies will be different for every hadron 
from the ones given by the dispersion relations, and in particular the rest energy of the neutral 
pion will not be zero (in the chiral limit). 

We can infer from our perturbative calculations, eq. (|2ip . that this renormalization is 
governed by one parameter, 02(59)1 with a further explicit dependence on a. It remains to 
be investigated whether this functional form is preserved at higher loop order and also non- 
perturbatively, and what will be the practical prescriptions that one has to infer from it. 
What is clear is that in numerical simulations the subtraction would depend on (3 and a in a 
different way from what the Callan-Symanzik renormalization group equations dictate in the 
scaling region. This derives from the further explicit dependence on the lattice spacing of the 
term multiplying C2. 

We end our conclusions with a discussion of a possible strategy to determine the renor- 
malized momentum non-perturbatively. If we denote the unsubtracted momentum by 
then the ir° correlation function at large Euclidean times will be 

^e^(vr (x)vr ot (0)) oc Ae"^ , x -» 00, (43) 

X 

where po is a function of the bare quark mass, m,Q, and the injected three-momentum p. The 
relation to the temporal component of the renormalized momentum, p^rem is then given by 

Po (f, m ) = po- ien (p, p c *;m ), (44) 

where denotes the value of the four-momentum for which the energies of the hadrons are 
restored to their physical values. Since the tt° is exactly massless in the chiral limit, the above 
relation can serve as a renormalization condition to fix the value of p c *. To this end, one has 
to evaluate po(p;mo) f° r several values of thq and p and determine p c * by implicitly solving 
the relation 

lim Po(p;m ) =Po- (45) 

m — >0 p=pfc r 

The energies in all hadronic channels can then be computed by setting the momenta to pjf. 
The technical difficulty associated with this strategy will be to control the large statistical 
fluctuations in the 7T° channel which arise from disconnected diagrams. 

One may wonder if other variants of minimally doubled fermions, such as the proposals 
by Karsten [9] and Wilczek [10] are easier to implement in simulations. We are currently 
investigating this issue, but leave a more thorough discussion to a future publication [25] . 
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